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Abstract 

We explore the capability of the IceCube/Deepcore array to discover signal neutri- 
nos resulting from the annihilations of Supersymmetric WIMPS that may be captured 
in the solar core. In this analysis, we use a previously generated set of ~ 70k model 
points in the 19-dimensional parameter space of the pMSSM which satisfy existing ex- 
perimental and theoretical constraints. Our calculations employ a realistic estimate of 
the IceCube/DeepCore effective area that has been modeled by the IceCube collabora- 
tion. We find that a large fraction of the pMSSM models are shown to have significant 
signal rates in the anticipated IceCube/DeepCore 1825 day dataset, including some 
prospects for an early discovery. Many models where the LSP only constitutes a small 
fraction of the total dark matter relic density are found to have observable rates. We 
investigate in detail the dependence of the signal neutrino fluxes on the LSP mass, 
weak eigenstate composition, annihilation products and thermal relic density, as well 
as on the spin-independent and spin-dependent scattering cross sections. Lastly, We 
compare the model coverage of IceCube/DeepCore to that obtainable in near- future 
direct detection experiments and to pMSSM searches at the 7 TeV LHC. 
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1 Introduction 



The nature of Dark Matter (DM), which makes up ~ 25% of the energy budget of the universe 
[H El E] is one of the greatest mysteries of modern physics. That DM is 'dark' implies that 
it is electrically neutral and thus has only been probed up to now through its gravitational 
interactions. Although there are many candidate particle physics scenarios which predict 
various kinds of DM j3], perhaps the most attractive possibility is that the bulk of the DM 
takes the form of a thermally produced, weakly interacting massive particle (WIMP) [3]. In 
most WIMP scenarios that also address the gauge hierarchy problem, the DM state is just 
the lightest, colorless and neutral member of an entire family of new particles; the lightest 
neutralino of R-Parity conserving supersymmetry (SUSY) [6] is arguably the most popular 
realization for DM in these kinds of models. 

Of course, at the end of the day, only detailed experimental measurements will be 
able to resolve the issue of DM's nature and data from multiple channels will be necessary. 
If DM is indeed a WIMP in a scenario such as described above, it could be discovered in 
the cascade decays of new, but heavier, colored states with TeV-scale masses which are 
produced with large cross sections at the LHC, with the DM appearing in the detectors 
as missing transverse energy (MET) [7J [5]. Similarly, the DM particle may scatter off of 
a nucleus, depositing energy in an underground detector and be directly observed via the 
resulting nuclear recoil [91 [10] • Finally, there are a number of ways that a WIMP may be 
indirectly observed through astrophysics experiments. For example, WIMP annihilation in 
the galactic center and/or halo, or in other galaxies, may be observed through the resulting 
annihilation products such as photons, positrons, antiprotons in satellite or ground-based 
observatories [HE]. A variant of this possibility is to observe the neutrinos that result from 
the annihilation of WIMPS that have been gravitationally captured in the core of the sun 
as it sweeps through the DM halo in the galaxy. This is the possibility that we will discuss 
below. 

In this paper we will explore the capability of the IceCube detector [T2], in conjunction 
with its densely instrumented sub-detector DeepCore [13J, which are installed in the ice near 
the South Pole, to detect neutrinos that would result from neutralino pair annihilation in 
the solar core. We derive our solar DM signals from models within the phenomenological 
Minimal Supersymmetric Standard Model (pMSSM), which is a 19-parameter version of 
the more general MSSM [H]. Specifically, we will examine the anticipated IceCube signal 
rates for a large number of points, ~ 70k, within this 19-dimensional space that lead to 
sparticle properties that are consistent with all of the existing experimental and observational 
constraints. We can then compare the capability of Ice Cube in exploring this 19-dimensional 
parameter space to that of other direct and indirect DM searches as well as to the ability of 
the LHC to find SUSY signatures within this same general framework. 

Why should one study such a large supersymmetric parameter space? General soft 
SUSY breaking within the MSSM leads to a scenario with over ~ 100 a priori free parameters 
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so that it is impossible to study in full detail. One approach in dealing with this problem is 
to consider specific, well-motivated SUSY breaking scenarios, such as mSUGRA [15J, GMSB 
[TU] or AMSB [T7] , all of which lead to a drastic reduction in the number of free parameters 
(to only ~ 3 — 5). Detailed studies of the resulting parameter space can then be easily 
performed. One problem with such studies is that they may bias us as to the nature of 
SUSY signals when performing DM (or collider) searches. An alternative to these top-down 
methods is to be less prejudicial and to instead follow a bottom- up approach such as that 
we have employed in a number of recent analyses [TSj [TUl [20], and will make further use 
of here. By imposing a set of theoretically and experimentally well-motivated constraints 
on the general MSSM, without making any reference to the specific mechanism of SUSY 
breaking, we are able to reduce the ~ 100 dimensional parameter space to one with 'only' 
19 parameters, i.e., the pMSSM, which is significantly more manageable. Such an approach 
has the advantage of being far more general than any specific SUSY breaking scenario and 
allows one to be in some sense agnostic about the SUSY mass spectrum, uncovering regions 
of parameter space that lead to distinct model characteristics and experimental signatures. 

There have been many previous studies of supersymmetric predictions for IceCube/ 
DeepCore (IC/DC), employing various strategies for simplifying the full MSSM parameter 
space. Detailed analyses of specific benchmark mSUGRA models were performed in [21]. 
In the work [22] predictions for points in 7- and 9- dimensional subspaces of the MSSM 
were investigated, while CMSSM predictions for IC/DC have been investigated in [39J and 
in a bayesian framework in [23]. Recent considerations of the IC/DC dark matter search 
outside of the context of the MSSM can be found in [21] [25]. In this work, we expand upon 
these previous studies and examine the capabilities of IC/DC to detect WIMP signatures 
from a large set of SUSY models (~ 70k) from the broad 19-dimensional parameter space 
of the pMSSM. We will find that a large fraction of these models have significant signal 
rates in the anticipated IC/DC 1825 day dataset. In some cases, the rate is large enough 
for early discovery. We find that observable rates are expected even for models where the 
LSP only constitutes a small fraction of the total dark matter relic density. In performing 
these analyses, we employ a realistic estimate of the IC/DC effective area as modeled by the 
IC/DC collaboration. In addition, we also examine in detail the dependence of the signal 
neutrino fluxes on the LSP mass, weak eigenstate composition, annihilation products and 
thermal relic density, as well as on the spin-independent and spin-dependent scattering cross 
sections. Coverage of the pMSSM model sample at IC/DC is compared to that obtainable 
in near-future direct detection experiments and to pMSSM searches at the 7 TeV LHC. 

In the next Section we will briefly discuss the techniques we employed in the gener- 
ation of the SUSY models corresponding to the ~ 70k pMSSM parameter space points we 
consider in this analysis and the various constraints that were imposed on their selection. In 
Section 3, we present our analysis method with the final results presented in Section 4. A 
discussion of these results and our conclusions will then follow in Section 5. 
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2 Generation of the pMSSM Model Set 



In this paper, we explore the sensitivy of the IceCube/Deep Core array in a broad region 
of Supersymmetric parameter space. This is similar in spirit to Ref. [19], which examined 
the indirect detection of dark matter in the pMSSM from its annihilation into electron- 
positron pairs. The set of pMSSM models that we investigate was generated previously in 
Ref. [13] and totals approximately ~ 70k points in the 19-dimensional pMSSM parameter 
space. Hereafter we refer to a point in this parameter space as a model. In this Section, we 
briefly review the procedure employed to generate this model sample. 

We study the 19-dimensional parameter space of the pMSSM [26]. This set of pa- 
rameters results from imposing the following minimal set of assumptions on the general 
R-Parity conserving MSSM: (i) the soft parameters are taken to be real so that there are no 
new CP-violating sources beyond those in the CKM matrix; (ii) Minimal Flavor Violation 
[27] is taken to be valid at the TeV scale; (Hi) the first two generations of sfermions with the 
same quantum numbers are taken to be degenerate and to have negligible Yukawa couplings 
and (iv) the Lighteset Supersymmetric Particle (LSP) is taken to be the lightest neutralino 
which is assumed to be a stable thermal WIMP. No assumptions about the physics at high 
scales or SUSY breaking mechanisms were employed. The first three of these conditions are 
applied in order to avoid issues associated with constraints from the flavor sector. These 
assumptions reduce the SUSY parameter space to 19 free soft-breaking parameters which are 
given by the three gaugino masses, Mj=i_ 3 , ten sfermion masses ^Q^Q^di^dsMMM^ 
the three A-terms associated with the third generation (A bit;T ), and the usual Higgs sector 
parameters /x, Ma and tan/3. 

To produce the set of pMSSM parameter points used in this paper, we performed 
numerical scans over the space formed by these 19 parameters. This required a selection of 
the parameter range intervals as well as an assumption about the nature of the scan prior 
for how the points were chosen within these intervals. These issues are described in detail 
in our previous works [TH [T51 [T9~] . Here, we simply note that two scans were performed: one 
employing a flat prior beginning with 10 7 points, and a second with a log prior employing 

2 x 10 6 points. The main distinctions between these two scans directly relevant to our 
analysis here are that (i) all SUSY mass parameters were restricted to be < 1 TeV for the 
flat prior case, while for the log sample the upper limit on mass parameters was raised to 

3 TeV, and (ii) the choice of the log prior generally leads to a more compressed sparticle 
spectrum than does the flat prior case. Note that the restriction on the upper limit for the 
mass parameters ensures relatively large production cross sections at the LHC for the case 
of the flat prior model sample. Most of the results presented below were obtained with the 
larger flat prior model set. 

Once these points were generated, we subjected them to a large set of theoretical and 
experimental constraints, and required consistency. This ensures that the model sets are valid 
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to study. We briefly review this set of restrictions here^j (z) The theoretical constraints we 
imposed are that the spectrum was required to be tachyon free, color and charge breaking 
minima must be avoided, a bounded Higgs potential must be obtained and electroweak 
symmetry breaking must be consistent, (ii) We employed a number of constraints from 
the flavor sector and precision electroweak data arising from the measurements of (g — 2) M , 
b — > S7, B — > Ti>, Bs — > meson-anti-meson mixing, the invisible width of the Z 

and Ap. (in) We required that the LSP contribution to the dark matter relic density not 
exceed the upper bound determined by WMAP. Note that we did not require the LSP to 
saturate the measured relic density; this leaves room for the existence of other dark matter 
candidates. In addition, limits from dark matter direct detection searches were also applied. 
(iv) The restrictions resulting from the numerous direct searches at LEP for both the SUSY 
particles themselves as well as the extended SUSY Higgs sector were imposed. Here, some 
care was necessary as some of these searches needed to be re-evaluated in detail due to SUSY 
model-dependent assumptions present in the analysis which we removed, (v ) Finally, the null 
results from the set of Tevatron SUSY sparticle and Higgs searches were imposed. The most 
restrictive data were found to be those from searches for stable charged particles [28J and 
those looking for an excess of multijet plus MET events |29j . We note that in the latter case, 
the search strategies were designed for kinematics expected in mSUGRA-inspired models. 
We thus simulated the search in some detail, at the level of fast Monte Carlo, for our full 
model sample. 

After this set of constraints was imposed, ~ 68.4fc models from the flat prior set 
survived this analysis chain, as well as a corresponding set of ~ 2.9k models from the log 
prior sample. This forms the set of models that we will consider in our following analysis of 
the IceCube/DeepCore capabilities to detect dark matter. 



3 Solar Neutrino Rate Calculations 

If the local DM halo is composed of WIMP dark matter with an empirically estimated local 
density po ~ 0.3 GeV/cm 3 we may be able to observe its presence as the result of WIMP 
scattering interactions in the sun. As the sun passes through the DM halo, WIMPs can 
scatter off of solar nuclei into orbits that are bound in the gravitational potential of the 
sun, eventually settling to the core of the sun after repeated scatterings [30]. Subsequent 
annihilations of captured WIMPs generically produce energetic particles, including neutrinos 
that can propagate out of the sun to terrestrial neutrino detectors. 

The instantaneous number of WIMPs in this captured population, N(t), can be mod- 
eled as 

^ = C c - C a N(t)\ (1) 

2 For full details, see Ref. Q3] 
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where C c is the capture rate and C a is proportional to the WIMP thermal annihilation cross 
section (av) (the constant of proportionality is the effective volume for WIMP scattering 
in the sun). The annihilation rate of captured WIMPs can be written as r a = C a N 2 . The 
contribution to the total capture rate C c from a shell of solar volume V at radius r is, 
schematically, 

dC - _ f j. tf( „, ^ W ai {v)-^v\ m(r)Pi(v,r). (2) 



J dvf(v,r)Y^ 



dV J J v ' ' ^ V y ' M y , 

j \ A. / 

The sum runs over the elements present in the sun. The number density of a specific 
target element is n.j(r) and the associated cross-section for WIMP-nuclei elastic scattering 
on target nuclei of this type is <jj. The probability of a scatter resulting in a bound orbit is 
Pi(v,r). The local sun frame WIMP velocity distribution is f(v,r) and p x is the local halo 
density of the scattering WIMP species^] Analytic expressions for these contributions can 
be found in [31]. For the sun, the dominant contribution to the capture rate is typically from 
spin- dependent elastic scattering of WIMPs off of hydrogen nuclei, (Tsd,p- Spin-independent 
elastic scattering off of heavier elements, though rare because of their low estimated solar 
abundances, may also provide important contributions due to the A 2 coherent enhancement 
of spin-independent scattering (oxygen, helium and neon may be most important in this 
regard [32]). 

Equation [T] can be easily solved; the annihilation rate of captured WIMPs at present 
time is given by [22J 

r Q ^^C a iV 2 (r ) = ^tanh 2 ^, (3) 

Z Z T e q 

where r eq = (C a C c )~^ is the time required for WIMPs to reach equilibrium and r = 
4.5 • 10 9 yr is the age of the sun. For r 3> r eq the solar WIMP population is in equilibrium, 
dN/dt = 0, and the annihilation rate approaches half the capture rate, 2r a /C c ~ 1. In 
this limit the annihilation rate has no explicit dependence on (av) and is instead completely 
determined by the capture rate, C c . The condition of equilibrium is commonly attained for 
the combinations of (av), asD, P and asi, P that are usually found for MSSM models, but is 
not completely general. This will be demonstrated in the more detailed discussion of the 
calculation of neutrino rates from our pMSSM models, to which we now turn. 



3.1 Solar Neutrino Signals from the pMSSM 

We now remark on some specific aspects of this calculation as they relate to our pMSSM 
model results and their interpretation. 

The magnitude and energy spectrum of the signal neutrino flux will depend on the 

1 p x may or may not equal the empirically estimated total DM density po ~ 0.3 GeV/cm 3 , as we describe 
in the next subsection LSP neutralinos in our pMSSM model set have p x < po, and the resulting neutrino 
signals must be appropriately scaled for each model. 
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properties of not only the LSP neutralino Xii but to some extent on substantially all of the 
masses and couplings needed to describe a particular pMSSM model. In order to calculate 
signal neutrino rates from each of our ~ 71k pMSSM models we rely on the computational 
package DarkSUSY 5.0.5 [22] • For each of our models we input SUSY Les Houches Accord 
files via SLHAlib [31] to DarkSUSY, which subsequently calculates annihilation and scatter- 
ing cross-sections. DarkSUSY uses the analytic formulae in [31] (as described in Ref. [22J) 
to derive the rate and spectra of signal neutrinos injected at the solar core by the solar Xi 
population. These injection spectra are then propagated out of the sun and to/through the 
earth to the detector via the package WimpSim [35]) which is embedded in DarkSUSY 5.0.5. 



It is important to note that, because we have employed the WMAP measurement 
of the DM relic density only as an upper bound on selecting pMSSM models, we must 
appropriately rescale the empirical estimate po = 0.3 GeV/cm 3 for the total local DM energy 
density by the factor 

R = (4) 
i 'WMAP 

where fiwMAP^ 2 = 0.1143 [36]. For definiteness we identify the set of pMSSM models which 
yield Q^oh 2 > 0.10 (R > 0.875) as models which saturate the WMAP bound and distinguish 
this subset of models in the figures that follow. LSPs in models with R substantially less 
than 1 are interpreted as comprising one component of a multi-component DM halo. 

As regards the local WIMP velocity distribution f(v,r), our calculations employ a 
standard Maxwellian distribution with dispersion v = 270 km/s (in the halo frame), boosted 
into the sun's frame (v = 220 km/s) and modified locally by the solar potential. It is difficult 
to quantify the uncertainty in the WIMP velocity distribution, but choosing a different form 
for the velocity distribution would likely yield at most an 0(1) change in the capture rate 
(see e.g. [22J). 

As mentioned in the previous section, solar WIMP capture-annihilation equilibrium 
has been seen to occur in most, but not all, of the previously studied MSSM models. We 
define capture-annihilation equilibrium for our models via the criterion: 

If > 0-9. (5) 



By this definition, 6.5% of our flat-prior and 14% of our log-prior pMSSM models are 
found to be out-of-equilibrium. Though out-of-equilibrium models typically lead to low solar 
neutrino rates we find a wide variety of predictions for out-of-equilibrium models in our set, 
and we will also distinguish this subset of models in the figures that describe our results. 

The numerical calculation of elastic scattering cross sections also deserves further 
comment. The calculation of x^-nuclei cross sections involves first the parton level calculation 
of Xi scattering off of the nucleonic constituents. Several different analytical treatments 
and associated numerical implementations of the parton-level scattering calculations are 
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available in DarkSUSY and the choice for this analysis is important because of the sheer 
phenomenological variety of the sparticle spectra in our pMSSM model set. We employ 
two non-default options: (i) to include neutralino and quark mass terms along with the 
squark mass terms in the denominators of propagators, and (ii) to include the DarkSUSY 
implementation of the extra twist-2 and box diagram contributions to the spin-independent 
cross section that were originally derived in [37] . The approximation to leave out these 
contributions is valid when squark masses are much greater than the quark and neutralino 
masses, which is not the case for generic models in our set. Overall, the default calculation 
was seen to diverge from more accurate calculations by more than an order of magnitude 
for ~ 20% of the models in our flat-prior set, and by up to five orders of magnitude for 
nearly pure bino-like LSPs which satisfied the WMAP relic density bound via large squark 
co-annihilations. Implementation of the non-default options addressed these issues. 

There are, of course, uncertainties that are also incurred in the estimates of the nuclear 
matrix elements that are necessary for going from the parton level cross-sections to the Xi~ 
nucleon cross sections (for recent accounts see [38]). These ambiguities have been shown 
to translate to an uncertainty in the capture rate of ~ 10% for dominantly spin-dependent 
contributions and by an 0(1) factor for dominantly spin-independent contributions [39 J. An 
additional ~ 20% error also applies to the latter case depending on the choice of solar models 
with differing metallicities [39J. In this work we use the following nuclear form factors (in 
the language of, e.g., [30]): 



fp 

■J u 


= 0.023, 


fP 

Jd 


= 0.033, 


fp 

J u 


= 0.26, 


rn 
J u 


= 0.018, 


rn 
Jd 


= 0.042, 


rn 
J a 


= 0.26, 


A£ 


= 0.842, 


A^ 


= -0.427, 


As 


= -0.085 



Other sources of uncertainty in the calculation include the modeling of the effects 
of charged and neutral current interactions as well as oscillations that neutrinos undergo in 
their propagation through the sun, to the earth and through the earth. We use the default 
DarkSUSY/WimpSim settings for solar composition, and neutrino propagation/oscillation 
parameters. The effects of varying these parameters on the resulting signal flux have been 
investigated in [35J, where it was found that significant uncertainties (~ factors of a few) 
may be present in the treatment of oscillations, particularly from our lack of knowledge of 
#13 and of the neutrino mass hierarchy. Accurate accounts of these errors will, of course, be 
important in determining the significance of any experimental result. However, as we are 
here focused on the SUSY model dependence of our predictions, we take into account the 
sources of error that we estimate to be the largest in this regard. 
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3.2 The IceCube/DeepCore Solar WIMP Search 



The IceCube neutrino telescope observes the Cerenkov light emitted from charged leptons 
as they transit through the Antarctic ice. The DeepCore extension of IceCube is a densely 
instrumented subsection of the larger IceCube detector, which is situated in the deepest ice 
near the center of the IceCube array [13]. This setup, using the larger IceCube detector as 
a veto for the embedded DeepCore detector, allows for a reduced neutrino energy threshold 
and hence better performance in searching for the neutrino signal resulting from WIMP 
annihilation. As neutrinos from WIMP annihilation pass through the ice and earth in and 
below the detector they produce detectable muons via charged current scattering which can 
then be identified and measured by the instrument. There are three backgrounds for the 
signal neutrinos produced by WIMP annihilation in the sun: (i) atmospheric muons, (ii) 
atmospheric muon neutrinos, and (Hi) neutrinos from cosmic ray interactions in the solar 
atmosphere. A large potential background from muons produced in cosmic ray interactions 
in the atmosphere is essentially eliminated by requiring candidate signal muons that come 
from the northern hemisphere, and hence from the conversion of muon-neutrinos that are 
capable of transit through a large portion of the earth. Because of this, the search that we 
describe here requires 1825 days (~ 5 yr) of useable data, but must be taken over a ~ 10 yr 
timescale. 

The number of observed signal events depends on the energy- and angular-dependent 
effective area of the detector. The effective area is determined by the IC/DC collaboration 
via monte carlo simulation. Here an initially generated population of neutrinos of various 
energies, incident on the earth from various directions, are transported through the earth, 
converted to muons, and transported to the detector, where, finally, a number of low-level 
data quality criteria are required in order to accept the event. In this way an energy- and 
angular-dependent efficiency is derived (i.e., taking the ratio of accepted events to generated 
neutrinos) which, when multiplied by the fiducial area of the detector, gives an effective area 
that should be convolved with raw signal neutrino fluxes^] in order to calculate the number 
of observed events. The effective area increases rapidly with increasing incident neutrino en- 
ergy as both the v^ — // conversion cross-section and the range of the produced muons grow 
approximately linearly with the neutrino energy; signal event rates are thus larger for stiffer 
signal neutrino spectra. The angular dependence of the effective area is fairly flat for neu- 
trinos incident from the northern hemisphere and it is dramatically suppressed for neutrinos 
incident from the southern hemisphere, where stringent angular cuts are necessary to reduce 
the atmospheric muon background. In performing the detector effective area convolution we 
integrate over the northern hemisphere only and over neutrino energies > 10 GeV. The raw 
neutrino fluxes given by the default DarkSUSY/WimpSim setup represent an average over 
the incident rates during the period spanning the spring and autumn equinoxes (i.e., the 
period during which the sun is below the horizon for a detector at the South Pole). 

3 In particular we use the "SMT8/SMT4" effective areas, derived as a preliminary estimate for the IceCube 
80-string+DeepCore configuration by the IC/DC collaboration [41] . 
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We define the detected neutrino rate, $^ (events/yr), as 

where d^ Ull}V/1 /dE are the raw differential neutrino/antineutrino fluxes calculated here using 
DarkSUSY/WimpSim, i£ is integrated over the neutrino energies > lOGeV, and A^' Vli (E) 
are the energy-dependent neutrino effective areas appropriate for neutrinos and antineutri- 
nos, respectively. 

This effective area does not include cuts which are needed to further reduce the 
atmospheric muon background and may lead to a reduction of the effective areas used here 
by a factor of ~ 5 [12] . To distinguish the annihilation signal from the isotropic atmospheric 
muon neutrino background, the IC/DC collaboration employs a directional search that looks 
for a statistical excess correlated with the solar angle. Ref. [12] estimates that the planned 
1825 day dataset will have a 90% exclusion sensitivity for a signal event rate of ~ 8 events/yr. 
Taking into account the factor of five signal reduction by implementing directional cuts, we 
estimate the 90% confidence exclusion limit for SMT8/SMT4 trigger level signal rates at 
$^ ~ 40 events/yr. From this we obtain a discovery threshold for trigger level signal 
events of $^ ~ 90 — 100 events/yr. In the work [23], it is estimated that 11.5 events/yr 
is the appropriate level for 90% exclusion by the Feldman- Cousins construction |43j . and 
that 31.6 events/yr is appropriate for discovery with S/\/B > 5. The range ~ 10 — 100 
events/yr thus constitutes a plausible estimate for exclusion/discovery criteria (where at 
the lower end one approaches the irreducible directional background of 0(10) trigger level 
events per year that are estimated from neutrinos produced in cosmic ray interactions in 
the solar atmosphere |44j |42j). More accurate criteria await a more detailed estimate from 
within the IceCube/DeepCore Collaboration. Throughout the rest of this paper we will use 
$^ = 40 events/yr as a criterion for exclusion, and when comparing to the discovery 
potential of LHC searches, we will use $^ = 100 events/yr as a conservative criterion for 
discovery. We note that the qualitative conclusions presented in the sections that follow 
are fairly insensitive to this choice, though the overall number of models that may be visible 
in our set varies considerably with the choice of discovery criteria. 

We note that the IceCube collaboration has also computed expected exclusion limits 
on WIMP-nucleon elastic scattering cross sections that may be attained in the 1825 day 
IC/DC search [15]. These limits are based on the conversion factors developed in [22J and 
assume that (i) the WIMP is in solar capture/annihilation equilibrium, (ii) the WIMPs make 
up the entirety of the local dark matter density and (Hi) the WIMPs annihilate exclusively 
into hard channels (i.e., W + W~ , Z°Z Q and r-pairs). We will compare this result with 
those obtained for the models in our set in Section 14.31 We note here that the authors 
of [22] discussed that differences in annihilation final state channels could generate about 
an order of magnitude variance in the predictions of models with fixed C c and m%o, and 
that, in addition, about one order of magnitude variance could be attributed to varying 
the LSP mass with fixed C c and annihilation final states (both effects are primarily due to 
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the highly energy-dependent detector effective area). This finding is largely echoed in our 
analysis. Since all of the other uncertainties that we have discussed are 0(1), we believe 
that our treatment of the detector effective area should allow us to take directly into account 
the dominant source of approximation, while allowing a more accurate interpretation of the 
SUSY model dependence of our results. 

4 Results from the pMSSM Model Set 

In this section we describe the detected neutrino flux rates that can be expected in Ice- 
Cube/DeepCore (IC/DC) for the models in our pMSSM set. We investigate the dependence 
of the resulting flux on relevant SUSY model parameters: m^o, asi tP , &sd, p , ^h 2 \^o and the 
annihilation rates into various SM final states. We compare the results obtained from models 
in our flat-prior set with those in our log-prior set. We combine estimates of background 
event rates with signal flux predictions to estimate which of our pMSSM models may be de- 
tected or ruled out by IC/DC in the 1825 day solar WIMP search data. Finally, we discuss 
the ability of IC/DC relative to current /planned direct detection experiments, as well as the 
LHC, to discover or constrain supersymmetric DM. 

4.1 Basic Results 

In Figures [Tjjs] we investigate the dependence of the detected neutrino flux, $^ (the result 
after convo luting the raw flux spectra with the detector effective area, Eqn. ([7])), on LSP 
mass, WIMP-proton elastic scattering cross-sections and LSP relic density. In each figure we 
display points for each of the pMSSM models in our flat-prior model set (grey) or log-prior 
model set (black). We highlight pMSSM models which are out of solar capture/annihilation 
equilibrium (as defined in Eqn. ([5])) in orange and pMSSM models that make up substantially 
all of DM, Vth 2 \ 

lsp ~ |wmap, in blue. In figures involving the elastic scattering cross 
sections, we use scaled cross sections, e.g., ctsd, p x R, where R was defined in Eqn. Q. Scaled 
cross-sections are appropriate because the limits placed on elastic scattering cross sections 
are proportional to the abundance of the WIMP and are usually quoted with the assumption 
that the scattering particle makes up all of the observed DM^j 

One very obvious feature in Figs. [3]j4] is that large elastic scattering cross-sections, 
especially the spin-dependent WIMP-proton cross-sections, are strongly correlated with large 
neutrino signal rates. This is to be expected as the normalization of the neutrino spectra 
from models in capture/annihilation equilibrium is determined by the capture rate, which is 
in turn determined by the elastic scattering cross-sections. We will turn to a more detailed 

4 In using this scaling, we are also assuming a canonical thermal cosmological history and that the DM 
distributions responsible for the signals are reasonably approximated by the large scale average abundance, 
i.e., that substructure in the DM distribution does not heavily affect the resulting signals. 
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Figure V. We display points representing models in our flat-prior pMSSM model set in 
the detected neutrino flux, <E>^ (see Eqn. Q), vs. LSP mass plane. Grey points represent 
all of the models in this set while orange points denote pMSSM models which are out of 
capture/annihilation equilibrium according to Eqn. ^ and blue points represent models 
whose LSPs form substantially all of dark matter, fl^,oh 2 > 0.10 (R ~ 1). 




Figure 2: We display points representing models in our flat-prior pMSSM model set in 
the detected neutrino flux, $^ (see Eqn. ([7])), vs. R (see Eqn. Q) plane. Grey points 
represent all of the models in this set, orange points denote pMSSM models which are out 
of capture/annihilation equilibrium according to Eqn. ^ and blue points represent models 
whose LSPs form substantially all of dark matter, Q^oh 2 > 0.10 (R ~ 1). R ~ 1 models are, 
of course, situated at the very right-hand edge this figure. 
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Figure 3: We display points representing models in our flat-prior pMSSM model set in 
the detected neutrino flux, $^ (see Eqn. Q), vs. asi, p xR plane. Grey points represent 
all of the models in this set, orange points denote pMSSM models which are out of cap- 
ture/annihilation equilibrium according to Eqn. ^ and blue points represent models whose 
LSPs form substantially all of dark matter, fi^o/i 2 > 0.10 (R ~ 1). We note that we include 
all models which satisfied the direct detection limits that were in place when these models 
were generated [T4] (wherein a factor of 4 error was allowed for nuclear form factor uncer- 
tainties). We will consider the effect of recent direct detection limits (including the very 
recent XENON100 limit \M) in Section |4~3 
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investigation of the SUSY model dependence of these results in the next section (Section 



From Figures [T]j4] it is difficult to tell how many pMSSM models would have high 
signal neutrino rates at IC/DC. Figure [5] includes a histogram of the detected neutrino flux 
rates in our flat-prior model set and we note here that approximately 83%, 48%, 8.6% and 
0.6% of these models are expected to have detected signal rates greater than 1, 10, 100, 
and 10 3 yr -1 , respectively. According to the background and signal significance estimates 
described in the previous section, these results imply that the IC/DC 1825 day dataset may 
be expected to exclude (discover) ~ 22% (8.6%) of the models in our flat-prior model set 
(the ~ 400 models with signal event rates > 10 3 yr -1 would likely be seen with much less 
than the planned 1825 days of data). Among the subset of pMSSM models whose LSPs make 
up substantially all of DM the coverage is somewhat better, with ~ 38% of such flat-prior 
models being excluded in the 1825 day IC/DC dataset. These results are displayed in Fig. 
[6] As expected, pMSSM models which are out of capture/annihilation equilibrium are much 
more difficult to exclude. Essentially none of the flat-prior out-of-equilibrium model subset 
are expected to be excluded by our criteria, as is demonstrated in Fig. |6j We also note that 
capture/annihilation equilibrium is not strongly correlated with relic density (this is quite 
insensitive to our definition of equilibrium) , see Fig. [2j 




Figure 4: We display points representing models in our flat-prior pMSSM model set in 
the detected neutrino flux, $^ (see Eqn. ([7])), vs. <jsd, p x R plane. Grey points represent 
all of the models in this set, orange points denote pMSSM models which are out of cap- 
ture/annihilation equilibrium according to Eqn. ^ and blue points represent models whose 
LSPs form substantially all of dark matter, Q^oh 2 > 0.10 (R ~ 1). We note that asD,p x R 
is the quantity which is most strongly correlated with high detected neutrino rates. 

In Figure [5] the detected neutrino rates for our log-prior pMSSM model set are com- 
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Figure 5: Comparison of the basic results for our flat-prior pMSSM model set with those 
for our log-prior pMSSM model set. Flat-prior models are represented by grey points while 
log-prior models are represented by black points. There are ~68.4k models in the flat-prior 
set and ~2.9k models in the log-prior set. We note the general similarity between these two 
model samples. 
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Figure 6: Using the approximate exclusion criterion described in Section ^2 we compare 
histograms of models that are expected to be excluded by IC/DC (red) against the full 
model set (grey) as functions of the LSP relic density and of the ratio 2T a /C c (recall Eqn. 
(|5])). As described in the text, exclusion performance is not strongly correlated with LSP 
relic density and we see that IC/DC may be able to discover LSPs with i? <C 1. In contrast 
to this we see that essentially none of the models that can be excluded are out-of-equilibrium. 
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pared with those for our flat-prior pMSSM model set. By the criterion described above we 
may expect to exclude about 25% of the models in our log-prior model set. We note that, 
besides the difference in statistics, the distributions shown in the histogram Fig. [5] are very 
similar to each other. While histograms like these are better interpreted as a description of 
the possibilities within our ~ 71 k pMSSM model set than as posterior distributions for the 
MSSM we note the encouraging robustness of these results with respect to a change in scan 
priors. 



From Figures [Tjj4] and Figure [6] it is interesting to note that there is not a strong 
correlation between the relic density of the LSP and the resulting neutrino signal. IC/DC is 
sensitive to pMSSM models with R > 10~ 2 at approximately the same level as those with 
R m 1. This is a reflection of the fact that, for models in capture/annihilation equilibrium, 
the solar WIMP signal scales with relic density]^ as $^ ~ {<Jsd, p x R) ~ 1, fc> r R ~ ^/( av )i 
similarly to the expected scaling in terrestrial direct detection experiments. This is in con- 
trast to the (av)R 2 ~ R scaling appropriate for the cosmic ray signals from DM annihilation 
in the galactic halo, galactic center (including, for example, IceCube searches for neutrinos 
from DM annihilating near the Galactic Center), or other astrophysical DM distributions. 
A more detailed discussion of the performance of IC/DC as it compares to/complements the 



LHC, direct detection and indirect detection experiments will follow in Section 4.3 



4.2 SUSY Model Dependence 

In order to examine the SUSY model dependence of these results, we first discuss the SUSY 
model dependence of the normalization of the signal neutrino spectra, and then we explore 
the dependence of the results on the shape of the signal neutrino spectra. 

For models in equilibrium, the normalization of the raw flux spectra is determined by 
the capture rate, C c . This relationship is illustrated in Figure [7j where we display our models 
as points in the detected neutrino flux vs. capture rate plane. As expected, we see that the 
detected signal fluxes are well correlated with the capture rates for models in equilibrium 
(the remaining scatter of the points reflects the variance in the shape of the signal spectra), 
while the correlation is not as strong for out-of-equilibrium models. The capture rate is, in 
turn, determined by a combination of asi,p, o~sd, P) ^ 2 |lsp an d rn^o. Figure [8] illustrates the 
correlation between the capture rate and elastic scattering cross-sections. 



In Figure [8] it is apparent that the correlation between C c and o~sd,p is stronger than 
that between C c and asi,p- This can be understood from Figures [9 10 as follows. In Figure 



9jwe observe the performance of the IC/DC search in the o~sd, p x R vs. asi, p xR plane. We 
display models that would be expected to be excluded by the IC/DC analysis (red points) and 
also those that we would not expect to be excluded (blue points) in this plane and see that 



5 This simple scaling is of course altered when the (av) effective at freeze-out is not the same as the 
effective (av) in present-day DM halos (for example when co-annihilations are significant [17]). Many of the 
models in our set display this property, though it does not significantly change the intuition of these remarks. 
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points that can be ruled out locate along approximate, but discernible, vertical and horizontal 
boundary lines on the figure. These lines roughly describe a reach of asD, P x R ~ 10~ 5 pb 
and <Jsi tP xR ~ 10~ 8 pb in this plane, and the ratio of these numbers implies that the IC/DC 
search is sensitive to spin-independent elastic scattering cross-sections that are roughly ~ 10 3 
times smaller than the thermal spin-dependent elastic scattering cross-sections to which it is 
sensitive (of course, terrestrial direct detection experiments currently probing asi, p xR have 
already surpassed this level of sensitivity). This ratio is essentially the ratio of the coefficients 
in Equation ([2]), which are, roughly, integrals over the solar volume of the abundances of 
various target nucle^J The question becomes: what are typical values of the ratio o- S d, p /o~si, p 
in our pMSSM model set? The distribution of values for this ratio, in both our flat- and 
log- prior model sets, is displayed in Figure 10 We find that ~ 66% of flat-prior and ~ 86% 
of log-prior models have (Jsd, p /o~si, p > 10 3 , so that for most of our models spin-dependent 
LSP-hydrogen scattering is the dominant mode of capture in the sun. 
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Figure 7: We display points representing models in our flat-prior pMSSM model set in the 
detected neutrino flux, vs. solar capture rate, C c , plane. Purple points represent models 
that are in capture/annihilation equilibrium according to Eqn. [5] and orange points denote 
pMSSM models which are out of capture/annihilation equilibrium. 

The model dependence that remains in Figures [7j(8j and the blurring of the horizontal 
and vertical "lines" in Figure [9j is due to the shape of the signal neutrino spectrum, and 
its effect on the conversion of raw flux spectra to detected flux spectra. This conversion is 
performed as discussed above, by convolving raw fluxes with detector effective areas that are 
sharply energy dependent (decreasing with decreasing energy). The shape of the spectra is 

6 The precise size of this ratio is, of course, dependent on the specific choice of solar composition model. 
Uncertainties due to the solar model were discussed in [39 and, while important in making accurate predic- 
tions, are not large enough to change the intuition of this discussion. 
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Figure 8: We display points representing models in our flat-prior pMSSM model set in the 
capture rate, C c , vs. elastic scattering, <tsi,p x R ot <tsd, p x R planes. Grey points represent 
all of the models in this set, orange points denote pMSSM models which are out of cap- 
ture/annihilation equilibrium according to Eqn. ^ and blue points represent models whose 
LSPs form substantially all of dark matter, fi^o/i 2 > 0.10 (R ~ 1). The strong correlation 
between C c and a S D, P x R explains the strong correlation observed between <E>^ and <jsd, p x R 
(Fig. [i]). The strong correlation between C c and usd,p x R (and relatively weak correlation 
between C c and &si,p x R) is explained by relative importance of spin-dependent scattering 
in solar WIMP capture and the fact that (Jsd,pI&si,p > 10 3 for most of the models in our set 
(cf. Figs. g|l0j. 



17 



3\ fl >40 Yr -1 Mods (Red) and $ V D <40 Yr" 1 Mods (Blue) 



-2 




-14- 

-16 -14 -12 -10 -8 -6 

Log 10 [tr S i, p *R (pb)] 

Figure 9: SUSY model dependence in the asD, P x R vs. a SIp xR plane. Flat-prior models 
that are estimated to be excluded in the IC/DC search are displayed as red points, those 
that are not expected to be excluded are displayed as blue points. 



obviously highly dependent on pMSSM model details, most directly via the DM annihilation 
rates into various SM final states. 



In order to study this dependence we focus on the ratio: 

(DetectedFlux) _ _ J A£(E)^ + A v £(E)^dE 

(RawFlux) -($ +$ ) + 



(8) 



In Figure 11 we display points in the (DetectedFlux)/(RawFlux) vs. LSP mass plane for each 
pMSSM model in our flat-prior model set (grey), and for models that annihilate with various 
levels of purity (described in the caption) into particular SM final states. This ratio varies by 
about two orders of magnitude over this model set, generally increasing with increasing LSP 
mass, as the resulting neutrino spectra are shifted in energy toward larger effective areas. 
There is about an order of magnitude spread in the ratio of (DetectedFlux) / (RawFlux) at 
any given LSP mass that is due to varying annihilation final states. As expected, the upper 
"edge" of the scatter is made up of pMSSM models annihilating dominantly into so-called 
"hard" channels, which yield relatively stiff neutrino spectra, i.e. W + W~, Z°Z° and r-pairs. 
The lower edge is populated by models annihilating dominantly to bb and thus resulting in 
relatively soft neutrino spectra^ The correlation between LSP eigenstate composition and 



7 The models with flux ratios that exceed the upper edge in the window m^o ~ 55 — 85 GeV are examples 
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Figure 10: We describe the distribution of spin-dependent and spin-independent WIMP- 
proton cross-sections in our flat- and log-prior model sets. We display points in the detected 
neutrino flux, vs. <JsD, P /ffsi,p plane. We also histogram the ratio <7sd, p / >cr si,p for both 
model sets. Flat-prior models are described by grey points/bars while log-prior models are 
described by black points/bars. We observe that most of the models in both sets have 
crSD,p/ffsi,p > 10 3 (the small subset of models for which asD, P < °~si, P are seen to have 
contributions to <Jsd,p from Z° and squark diagrams which largely cancel). 
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this ratio, (DetectedFlux)/(RawFlux), is not as strong as that between annihilation final 



states and this ratio (cf. Figs. 11 and 12) as the connection between the two is muddied 
somewhat by the details of the superpartner mass spectrum. Despite this, some amount of 
correlation is still evident. Nearly pure higgsino and wino LSPs typically annihilate with 
large rates into the W + W~ and Z°Z° final states, thus tending to populate the upper edge 
in the figures. Nearly pure bino LSPs annihilate dominantly through helicity suppressed 
sfermion exchange graphs so that (for m^o < m t ) annihilations of bino LSPs proceed to 
a mixture of bb and r-pair final states, and the resulting (DetectedFlux) / (RawFlux) ratios 
interpolate between the pure r and pure b results. 



In Figures |13|fl4| we show the detected flux histograms for various annihilation final 
states and LSP eigenstate compositions. One particularly interesting result is that the 
detected fluxes for mixed neutralinos are found to be typically quite high. This can be seen as 
the consequence of several factors. The most important factor in attaining high fluxes is the 
combination (TsD,pXR- The spin-dependent cross-section itself is due to Z° or squark mediated 
graphs. In any LSP eigenstate scenario the squark mediated contributions to scattering 
exist and are largely dependent on the scanned squark masses. However, potentialljj^] larger 
cross-sections exist in cases where the LSP has significant higgsino fraction, and thus can 
couple strongly to the Z°. As a competing effect, a large higgsino content opens several new 
annihilation channels and has been shown in simplified models to give LSP relic densities 
that fall significantly short of the WMAP measured value (for LSPs below ~ 1 TeV [IE])- 
The extent to which lower relic density implies a lower overall value of <j,sd,p x R depends 
upon how closely the present day elastic scattering cross-sections are related to the freeze- 
out annihilation cross section. Annihilation sub-processes such as XiXi ~^ K H, Z° — > qq 
decrease the relic density while simultaneously (potentially) increasing SI and SD scattering 
rates. In contrast, annihilations to non-quark final states only deplete the LSP relic density 
(resonant funnels are also far more efficient in s-channel annihilations than in t-channel elastic 
scattering). Finally, as noted above, higgsino-like LSPs may have significant annihilation 
rates into the W + W~ and Z°Z° final state channels, i.e., those which see the largest detector 
effective area. LSPs with sizeable wino content share this latter virtue (annihilation to hard 
channels), but because of the lack of couplings to h, Z°, as well as small xt~ Xi mass splittings 
which lead to efficient annihilation, they do not usually have high values of asn,p><R- A look 
into the subset of mixed LSPs in our flat-prior model set shows that essentially all models 
have a significant higgsino fraction (97% of the mixed models have LSPs that are > 20% 
higgsino), so that many of these models yield observable solar neutrino rates. 



of LSPs whose annihilations are largely to the jZ° final state, resulting in a neutrino spectrum from the 
Z° decay that is even stiffer than the spectra that would result from annihilation into the Z°Z° final state 
(these models nevertheless end up with small signal rates because of generally low relic density). 

8 In principle the Z° and squark contributions, with associated signed nuclear form factors, may add 
constructively or cancel against each other, depending on the details of the model spectrum and assumed 
values for the nuclear form factors. 



20 



All FLAT (Gray) 



(<crv>),/(<crv>)>0.85 




oil 
C 



-7 

y 



aa 

3 -10 



L 2 o 



50 100 200 
LSP Mass (GeV) 

(<crv>) r /(«jv>)>0.70 




50 100 200 
LSP Mass (GeV) 

(< trv > ) iy/(<cr v > ) > .99 




50 100 200 
LSP Mass (GeV) 



500 




-7 



T -8 



3D 

c 



10 





-7r 


e 


-8 




-9 


a 

c 




ol) 




5 


-10 




-4 



50 100 200 
LSP Mass (GeV) 

(<crv > V(<<xv>)>0 .93 





















■ 4 







^0 50 100 200 

LSP Mass (GeV) 

(<trv>) z /(<crv >)>0 .42 



500 




50 100 200 
LSP Mass (GeV) 



500 



Figure 11: We display points representing our flat prior models in the 
(DetectedFlux)/(RawFlux) (Eqn. rt8|) vs. LSP mass plane. The full flat-prior model set is 
displayed as grey points and models whose annihilations occur predominantly through a 
given final state channel are overlaid in other colors. Models with [av) t i/ (cry) > 0.85 (cyan), 
with (av) Tf / (av) > 0.70 (green), with (av) b i/(av ) > 0.93 (red), with (av )w+w- / {w) > 0.99 
(blue) and with (av)zoz°/(&v) > 0.42 (magenta) are shown. Purities are chosen to obtain 
subsets of models of similar size (in our model set there is a maximum purity for annihi- 
lations into the Z°Z° final state, (av)zoz /(w) < 0.445, as the higgsino-like LSPs which 
annihilate well to Z°Z° via x\ exchange also annihilate well to W + W~ via xt exchange). 
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Figure 12: We display points representing our flat prior models in the 
(DetectedFlux)/(RawFlux) (Eqn. (J8|) vs. LSP mass plane. The full flat-prior model 
set is displayed as grey points and models categorized according to LSP eigenstate compo- 
sition are overlaid in other colors. By convention our LSPs are described in terms of their 
neutralino mass matrix entries as: Xi = %uB + Z 12 W 3 + + Zi^H®. Higgsino models 

are defined as having (I.Z13I 2 + |Z 14 | 2 ) > 0.99 and are displayed here in green. Wino models 
are defined as having | 1 2 > 0.99 and are displayed here in blue. Bino models are defined 
as having \Zn\ 2 > 0.99 and are displayed here in red. Mixed models are defined as having 
|Zn| 2 , \Zi 2 \ 2 and (|Zi3| 2 + |Zi 4 | 2 ) all < 0.8 and are displayed here in magenta. 
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Figure 13: Detected neutrino flux histograms illustrating the dependence of $^ on final 
state annihilation channel. The full flat-prior model set is described in the grey histogram 
while subsets of these models whose annihilations occur predominantly through a given final 
state channel are displayed in other colors. Models with (av) t i/ (av) > 0.85 (cyan), with 
(av) Tf / (av) > 0.70 (green), with (av) h i/ (av) > 0.93 (red), with (av) w + w - / (av) > 0.99 
(blue) and with (av) z o z o / (av ) > 0.42 (magenta) are shown. Purities are chosen to obtain 
subsets of models of similar size (in our model set there is a maximum purity for annihilations 
into the Z°Z° final state, (av) z®z° / (crv) < 0.445, as the higgsino-like LSPs which annihilate 
well to Z Q Z° via x® exchange also annihilate well to W + W~ via xt exchange). 
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Figure 14: Detected neutrino flux histograms illustrating the dependence of $^ on LSP 
eigenstate composition. By convention our LSPs are described in terms of their neutralino 
mass matrix entries as: Xi = ZnB + Zi 2 W 3 + Zi 3 H® + ZnH®. Higgsino models are defined 
as having (|^i3| 2 + |-Zi 4 | 2 ) > 0.99 and are displayed here in green. Wino models are defined 
as having \Zi 2 \ 2 > 0.99 and are displayed here in blue. Bino models are defined as having 
\Zn\ 2 > 0.99 and are displayed here in red. Mixed models are defined as having \Zu\ 2 , \Zi 2 \ 2 
and (|^i3| 2 + \Zu\ 2 ) all < 0.8 and are displayed here in magenta. 
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4.3 IC/DC Comparison with Other Experiments 



We now compare the ability of IC/DC to constrain or discover the pMSSM models in our set 
with that of direct detection experiments, searches in other indirect detection experiments, 
and at the 7TeV LHC. In order to make this comparison, we focus on constraints in the 
asi, p xR vs. LSP mass, <tsd, p xR vs. LSP mass and asi, P xR vs. cjsd, p xR planes. In the 
following figures, pMSSM models which are expected to be excluded by the IC/DC 1825 day 
dataset {i.e., those with detected neutrino fluxes $^ > 40 yr _1 ) are displayed as red points 
over the grey points which represent the entire flat-prior model set. pMSSM models which 
cannot be excluded are displayed in blue. 

For comparison to current and planned direct detection experiments we combine the 
neutrino fluxes calculated in this work with quoted limits on the thermal spin-independent 
cross-section, a SI)P xR } from CDMS \4S\ and XENON100 j5U], as well as the recent XENON100 
limit [46] . We also show projected limits on <r S i p xR from LUX (30,000 kg/days) [51], Super- 



CDMS at SNOLAB [52] and COUPP (60kg and 500kg lyr searches [53]) in Figure (T5J Here, 
we recall that in generating our pMSSM model set [14], we employed the direct detection lim- 
its that were applicable at the time as a constraint on viable models (while allowing a factor 
of 4 error to account for the uncertainties associated with computing the spin-independent 
LSP-nucleonic scattering cross-section). The recent XENON100 result [46] represents the 
current most stringent limit on a,si, P xR and certainly rules out some subset of our pMSSM 
model set. While we have not dropped these models from the present discussion, we will 
consider the impact of this result in what follows. For the thermal spin-dependent scattering 
cross-section, &sd, p x R, we show the quoted limits from AMANDA (2001-2003 data) [43] 
and the IceCube-22 string (no DeepCore) 2007 data [54] , we also display projected limits on 
&sd, p xR from COUPP (4kg [55], 60kg and 500kg [53] lyr searches) and IceCube/DeepCore 



56] in Figure 16 Recall that limits from Amanda and IceCube are placed under the as- 



sumption of annihilation dominantly to either hard or soft final state channels. 

While there is substantial overlap between excluded and non-excluded models in both 
of the quantities asi, p xR and <jsd, p xR, we observe that essentially all of pMSSM models with 
asi, p xR > 10~ 7 pb or ctsd, p xR > 10~ 4 pb, and a majority of models with asi lP xR > 10~ 8 pb 
or a S D, P xR > 10~ 5 pb, would be excluded by the IC/DC search. The recent XENON100 



limit (the black-dashed curve in Fig. 15) is ~ 10~ 8 pb over most of the LSP masses in our 
set. This value excludes ~ 16% of the models in our flat-prior set, although we note that 
~ 10~ 8 pb lies on the steep portion of the <Tsi, P x R distribution for these models so that, if 
nuclear form factors were to conspire so that the appropriate limit were looser or tighter by 
a factor of 2, ~ 10% or ~ 25% of our flat-prior models would be excluded, respectively. If we 
assume, as a rough estimate, that all models with asi, P xR > 10 -9 pb will be excluded in the 
near-future ton-scale spin-independent scattering direct searches, we estimate that, in our 
flat-prior model set, ~18% of the models will be excluded by both the IC/DC search and the 
spin- independent direct detection (SI DD) searches, while ~4% will be excluded by IC/DC 
and not by SI DD searches, ~31% will be excluded by SI DD searches but not by IC/DC 
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and ~47% of these models would not be excluded by either IC/DC or SI DD searches. 

Recall that the previously estimated IC/DC exclusion limit (the black line in Figure 



16, based on the derivation presented in in |22j) assumes (i) annihilations exclusively to hard 
channels, (ii) an LSP in capture/annihilation equilibrium and (Hi) an LSP that makes up 
all of DM (R = 1). One can see that this line agrees fairly well with the putative upper line 
described by the set of our models which IC/DC cannot exclude (the blue points in Figure 
16); only a small number of such models have <JsD,pXR above this line. Upon examination 



these are models which annihilate dominantly to soft channels (violating assumption (i) 
above) and, for non-excluded models well above the line, a significant amount of squark 
coannihilation is present so that the scattering cross sections are relatively large (due to the 
light squark) but still difficult to constrain as R < 1 (violating assumption (ii) above). As 
expected there are a large number of models that may be excluded with &sD,pXR ranging 
much below this line: these are models for which C c is dominated by spin-independent 
scattering. A large fraction of these models would thus also be expected to be excluded by 
SI DD experiments. 

We note that the IC/DC and COUPP 60kg searches seem to do a similar job of 
probing both o~sd,p x R a> n d o~si, p xR, although it may not be obvious from the information 



provided in Figs. 15 16 alone. From Figure 16 there are apparently quite a few models with 
&sd, p x R far below the sensitivity expected by COUPP that may indeed be excluded by 
IC/DC. However, these are necessarily models for which C c receives a sizable contribution 
from asi,pXR as discussed above. We would expect that asi lP xR <; 1CT 9 pb for such models, 
roughly within the asi, p xR sensitivity expected of the COUPP 60kg search. We thus expect 
the two experiments to provide highly complementary cross-checks on any WIMP discovery 
or limits in the o~sd, p xR vs. a$L p xR plane. 

For comparison to potential discoveries at the LHC we employ the expected results 
of an ATLAS 4j01 (i.e., 4-jet + 0-lepton + missing transverse energy) search (y/s = 7TeV, 
1 fb _1 and 50% systematic uncertainty are assumed) that were found for this same set 
of pMSSM models in the work [57j^J Given that models characterized as "passing" or 
"failing" the ATLAS analysis was based on a discovery criterion, rather than an exclusion 
criterion, we compare these results to our (conservative) criterion for discovery in IC/DC, 



$^ = 100 events/yr (see Section 3.2). Using this criterion, we note that ~ 8.6% of the 
models in our flat-prior pMSSM set are expected to be discoverable in the IC/DC solar 
WIMP search. 



Figure 17 demonstrates the combined ability of the IC/DC solar WIMP search and 
the ATLAS 4j01 search to discover models in our flat-prior pMSSM model set. In this figure 
we display models that would be expected to be seen by both searches, as well as models 
that would be discovered at IC/DC and missed in the 4j01 search, models that would be 
missed at IC/DC and seen in the 4j01 search and finally models expected to be missed by 



9 Results for other search channels are similar. We use the 4j01 search channel here as it was found to be 
the most effective discovery channel for the pMSSM [57j . 
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Figure 15: Comparison of IC/DC and spin-independent direct detection searches. We display 
all points in the flat-prior model set in grey, models that are estimated to be excluded by 
the IC/DC solar WIMP search in red and those which are not excluded in blue. Current 
experimental limits from the CDMS and XENON100 [50J [46J collaborations are displayed 
as red (CDMS 2010), black-solid (XENON100 2010) and black-dashed (XENON100 2011), 
respectively. Near-future projected experimental limits from LUX [51J, SuperCDMS [52J, 
COUPP 60kg and COUPP 500kg [53] are displayed as brown, magenta, orange, green-dashed 
and green-dotted lines, respectively. 
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Figure 16: Comparison of IC/DC and spin-dependent direct detection searches. We display 
all points in the flat-prior model set in grey, models that are estimated to be excluded by the 
IC/DC solar WIMP search in red and those which are not estimated to be excluded in blue. 
Current experimental limits from the AMANDA [15] and IceCube-22 [51] collaborations 
are displayed as orange and magenta lines, respectively (with the assumption of soft or 
hard channel annihilations represented by dashed or solid lines, respectively). Near-future 
projected experimental limits from the COUPP [55J[53] 4kg, 60kg and 500kg searches in 
green- solid, dashed and dotted lines, respectively. The IceCube/DeepCore limit estimated in 
[56] (assuming hard channel annihilations, DM which is in capture/ annihilation equilibrium 
and DM which has R = 1) is displayed as a black line. 
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both searches in the <jsd, p xR vs. asi, p xR plane. Of course, in the time that will be required 
for IC/DC to perform the search described in this paper, the LHC will likely accumulate 
much more than 1 fb _1 of data. However, the expected coverage of our model set by LHC 
searches using , e.g., 10 fb -1 of data, is > 90% [57] (the vast majority of our pMSSM models 
would be seen at the LHC) and it makes little sense to ask which models are seen or unseen 
by combinations of IC/DC and the LHC in this scenario. The point of Figure 17 is not to 
compare coverage of the two classes of experiment at some definite time in the future, but to 
observe that the shape of these regions are determined primarily by the IC/DC search (cf. 
Fig. [9]). This is illustrative of the fact that the solar WIMP signal is very strongly correlated 
with the elastic scattering cross-sections while the rate of 4j01 events in ATLAS is, of course, 
only very indirectly so. 

As a final remark, we discuss the relative performance of searches sensitive to the DM 
annihilation cross-section, (av); a detailed study of such signals from DM in this pMSSM 
model set was investigated in [19J. As mentioned before, in contrast to searches in which 
signals are largely determined by elastic scattering cross sections (i.e., in direct detection 
experiments and the IC/DC solar WIMP search), indirect searches for DM annihilation 
will probe signals ~(crt')i? 2 . In lieu of non-standard cosmological scenarios [58] [59], which 
may drastically alter the connection between the annihilation cross-section and the relic 
density, or special relationships in the SUSY mass spectrum that allow for co-annihilations 
or resonant annihilations, we expect the scaling R ~ l/(av), so that (av)R 2 ~ R while 
o~sd,p x R ~ 1- Thus indirect searches for DM annihilation that probe (av)R 2 , such as the 
PAMELA /FERMI /AMS-02 cosmic-ray antimatter measurements, FERMI /MAGIC /HESS 
7-ray measurements and IC/DC observations of neutrinos from the Galactic Center, are 
much more sensitive to the LSP relic density than the other classes of experiments. Such 
experiments will typically have a much harder time discovering LSPs with R <C 1. We 
note that the largest annihilation cross sections in our flat-prior set are (o~v)R 2 ~ 6 x 
10 _26 cm 3 s _1 , and only about a quarter of the models in this set have {av)R 2 > 10~ 27 cm 3 s _1 . 
While astrophysical indirect detection limits typically come with much larger uncertainties 
than other classes of experiments, due to the difficulty in estimating the strength of the 
annihilation source (e.g., in estimating DM halo profiles and substructure [60]), a sensitivity 
at the level of (o~v)R 2 > 10~ 27 cm 3 s -1 or better, for canonical choices of profile, etc., will 
probably be necessary in order to be sensitive to a large fraction of models in our set. 



5 Discussion 



In this paper we have investigated the ability of the upcoming 1825 day IceCube/DeepCore 
solar WIMP search to discover /constrain SUSY WIMPs. In this aim we have employed the 
large set of ~ 71k phenomenologically-viable pMSSM SUSY models that were generated and 
described in the work [14]. We have discussed the basic calculation of neutrino telescope 
signals from captured WIMPs annihilating in the solar core and the details of our analysis, 
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Figure 17: Comparison of IC/DC and the ATLAS 4j01 (y/s = 7TeV, 1 fb" 1 and 50% sys- 
tematic uncertainty assumed) searches. Here "PASS" and "FAIL" denote discovered or not 
discovered, respectively. We display all points in the flat-prior model set in grey, models that 
are estimated to be discoverable by both the IC/DC solar WIMP search and the ATLAS 
search are displayed in pink, those expected to be seen by IC/DC but missed in the ATLAS 
search in green, those expected to be missed by IC/DC but seen in the ATLAS search in 
purple and those which are estimated to be unobservable by either search are displayed in 
black. 
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which relied heavily on the use of the computational package DarkSUSY [33] . 

We have discussed sources of uncertainty that affect the capture process. Many works 
have previously elucidated important sources of experimental uncertainty [38] [39] (i.e., from 
our imprecise knowledge of nuclear matrix elements, solar composition and of the details of 
the neutrino sector) in the estimation of the capture rate and subsequent neutrino signal. 
Although we made no attempt to account for such uncertainties here, we have found that, 
given the diversity of our pMSSM model spectra, much larger errors (orders of magnitude) 
in the resulting rates would result from a poor choice of parton-level scattering amplitude 
calculations and/or failure to take into account the IC/DC detector effective area. Both of 
these have been carefully considered in this work. 

As a basic result of our analysis, we find that a large fraction of our pMSSM models 
are expected to have significant signal rates in this search. We find that LSPs with a 
wide range of masses can be excluded and, somewhat surprisingly, the IC/DC search reach 
extends to many of the models for which the LSP forms only a small fraction of the total 
DM abundance. A study of solar capture/annihilation equilibrium confirms the expected 
result, that essentially none of the out-of-equilibrium LSPs in our set are are expected to be 
excluded. We have compared the results from subsets of pMSSM models that were generated 
with flat- and log-prior scanning over parameters and have found them to be quite similar. 

We have described the SUSY model dependence that is seen in the determination of 
both the shape and normalization of the resulting signal neutrino spectra. Since most all 
of the models that may be excluded by IC/DC are in capture/annihilation equilibrium the 
normalization of the spectra is essentially determined by the capture rate and since most 
of our models have ctsd,p/&si,p > 10 3 and the sun is largely composed of hydrogen nuclei 
targets, the normalization is mostly determined by <jsd,p x R- The shape of the spectra are 
of course most closely tied to the annihilation final state channels and we find that, for a 
given LSP mass, the final state channels producing the hardest and softest neutrino spectra 
see about an order of magnitude difference in the effective area of the detector. 

Given the relative importance of the SD cross-section, relic density and final state 
annihilation channel we find semi-predictive differences between the results expected from 
classes of models with bino, wino, higgsino and mixed LSP eigenstates. We observe that 
nearly pure wino or higgsino LSP models typically have low (though not always, unobserv- 
able) rates due to their generally very low LSP relic densities. The bino LSP models have 
the widest range of predictions as both their relic density and elastic scattering cross sections 
are largely determined by the scanned sfermion masses. We noted that mixed LSP models 
typically predict very large rates and understand this as arising from their LSPs almost al- 
ways having a significant higgsino fraction. Such models can thus attain large SD scattering 
cross sections and large annihilation rates into hard channel final states, without the very 
low relic densities that would result from being very purely higgsino. 

We conclude with our expectation that the IceCube/DeepCore search may play an in- 
tegral role in the experimental confirmation or expulsion of Supersymmetric neutralino dark 
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matter. In the event of production and measurement of new invisible states at the LHC there 
can be no guarantee of their cosmological stability or relic abundance until complementary 
observations can be made in direct or indirect detection experiments. While terrestrial direct 
detection experiments are already poised to probe deep into the space of scattering cross sec- 
tions populated by our models, the IC/DC search offers competitive sensitivity to the same 
scattering cross sections on a similar time scale and with somewhat orthogonal systematics. 
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